Charge and orbital order in RNiOs (R=Nd,Y) by LSDA+U method 
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Electronic structures of NdNiOa and YNiOa are calculated by using LSDA+U method with rota- 
tional invariance. The Jahn- Teller distortion is not allowed under the observed magnetic ordering. 
No orbital order on Ni sites cannot be observed in the calculation of both systems but different 
types of charge ordering. In a small distorting system NdNiOs, all Ni ions are trivalent and oxygen 
sites have a particular ordering, both in charge and orbital. In a large distorting system YNiOs, the 
charge disproportionation occurs, 2Ni^+ Ni^^+Ni*"*". Therefore, the charge ordering stabilizes 
the asymmetry of the arrangement of Ni magnetic moments in both systems. 
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Perovskite transition metal oxides are now of high in- 
terest because of the large variety and the possible con- 
trollability of their physical properties. Perovskite nicke- 
lates RNiOs (R = an trivalent rare earth or Y ions) can 
be classified into three different categories, according to 
their tolerance factor {t = dn-o / V^dm-o) ■ One is 
those whose tolerance factor is much smaller [t < 0.91) 
and has a larger distortion, such as LuNiOa or YNiOa. 
They are antiferromagnetic insulators at low tempera- 
tures. This class of RNiOs undergoes transition to para- 
magnetic insulator at the Neel temperature Tn and also 
has another phase transition from paramagnetic insula- 
tor to paramagnetic metal (M-I transition) at very high 
transition temperature. Second is those whose tolerance 
factor is intermediate (0.91 < t < 0.93), such as NdNiOa 
or PrNiOa. These materials are also antiferromagnetic 
insulators at low temperatures and above Tn they are 
paramagnetic metals. § Third is LaNiOa where t = 0.94 
and is paramagnetic metal. 

Low temperature phase of above mentioned first and 
second classes of RNiOa have the unique magnetic struc- 
tures. Its magnetic diffraction peak is characterized 
by the propagation vector k = (1/2, 0, 1/2) and the mag- 
netic unit cell is identified as the 2x1x2 supercell of 
the crystallographic lattice. [0^] The magnetic order is 
specified to be alternating ferromagnetic (FM) and anti- 

FM AFM 

ferromagnetic (AFM) couphngs •••NiT— Q — Ni^— O 

FM AFM FM . 

-Ni^- O -Ni^- O -Nil- O -Nil ••• along three 
crystallographic directions. 

On the contrary to the asymmetric arrangement of 
magnetic bonds (FM/AFM) around each nickel site, Ni- 
O distances around each Ni ion are almost equal to each 
other. For instance, the difference between the longest 
and the shortest Ni-0 bonds in one NiOg octahedron 
is at most 2% in RNiOs. ||l|,|,| The Jahn-Teller distor- 
tion is absent in RNiOs, even in the largely distorted 
system YNiOa, whose tolerance factor is 0.88. Exper- 
imentally there is only one crystallographic site for Ni 
ions in NdNiOa and two different sites for Ni ions in 
YNiOs. Therefore, one can expect trivalent ions Ni'^+ 



in NdNiOa . On the contrary, Ni ions in YNiOs are diva- 
lent and quadrivalent, i.e. the charge disproportionation 
2Ni3+ m^++m^+ occurs. 

In this letter, we study the electronic properties of 
the AFM insulating phase of two perovskite nickelates, 
NdNiOa and YNiOa, which are typical ones of above 
mentioned two respective classes. We use the LSDA+U 
method with the rotational invariance in conjunction 
with the LMTO-ASA method. § The LSDA+U method 
counts the electron-electron interaction between local- 
ized orbitals by the Hartree-Fock type interaction term. 
The ionic positions and lattice parameters used in the 
present calculations are imported from the diffraction 
experiments. |jl|,Q The unit cell here is the 2x1x2 
crystallographic supercell containing 16 Nd or Y ions, 
16 Ni ions, and 48 O ions and, in addition, 32 empty 
spheres, totally 112 atomic spheres. Sixty- four k-points 
in the Brillouin zone are sampled in the calculation of 
the density of states. NdNiOa is antiferromagnetic up 
to Tn ~ 200K, but the Nd spin moment vanishes at 
about 30K. |^ Therefore, three 4f electrons of Nd ion are 
counted in the frozen core. We also calculated the elec- 
tronic structure of NdNiOs with 4f electrons in valence 
states, and the results show no significant difference. The 
Coulomb and exchange parameters U and J of Ni ions 
are fixed to be 7.0 eV and 0.88 eV, respectively, through 
all the LSDA+U calculations. These values are consis- 
tent with photocmission experiments and the results of 
the LSDA calculations. |^,|| Calculated results are not 
changed much over a large range of values of U and J. 

The crystallographic space groups of NdNiOa and 
YNiOs are orthorhombic Pbnm and monoclinic P2i/n, 
respectively. The group theoretical analysis shows that, 
using the projection operator method, the Jahn-Teller 
distortion of one NiOe octahedron can not be trans- 
fered over the whole enlarged 2x1x2 supercell with 
k = (1/2,0,1/2). Therefore, the Jahn-Teller distor- 
tion actually cannot appear in the observed spin ordered 
state. 

NdNiOs: There are two possible spin configurations 
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in NdNiOs satisfying the observed transfer vectors k = 
(1/2,0,1/2). One is that the spin magnetic moments 
aUgn ferromagnetically on the (101) plane and the planes 
are stacked with a doubled period as TTiiTTii ? 
whose magnetic space group is monoclinic Pba. ||ll| This 
spin configuration is assumed in the present letter. An- 
other possible spin configuration with ttii spin order 
could be the doubled period checkerboard stacking along 
crystallographic y-axis, and the magnetic space group is 
orthorhombic Camc2i. |ll] The state of Camc2i shows 
orbital ordering of e^ states on Ni sites, but the resultant 
calculated total energy is higher by 0.18 eV per 2x1x2 
cell than that of Pf,a. Because there is no possibility of 
subsidiary Jahn- Teller distortion to lower the total en- 
ergy, CaTnc2i could not be the symmetry to be consid- 
ered. The total energy is also calculated in a fictitious 
antiferromagnetic state of Pb'nm' whose magnetic unit 
cell is identical to the crystallographic one and is higher 
by 2.0 eV per 1x1x1 cell. 



FIG. 1. Projected density of states of trivalent nickel Ni 
in NdNiOs. The energy zeroth is set at the top of the oc- 
cupied states Ep. The local z axis directs from Ni to one of 
oxygen atoms with the shortest Ni-0 distance, and the local 
X and y axes to the appropriate O atoms. 

Figure |l| shows the projected density of states on the 
Ni 3d orbitals in NdNiOs. The system is insulator with a 
gap Eg = 0.11 eV. Each nickel ion has the local magnetic 
moment ±1.1/1^ within the atomic sphere of a radius 
2.51 au, contrary to the observed value of ±0.9fiB- It 
should be noticed that the local magnetic moment cannot 
be uniquely defined. Furthermore, NdNiOs is not a sim- 
ple antiferromagnetic insulator but a dynamical effect is 
essential, which may be an origin of this discrepancy. |^] 
There is no distinctive variation in partial spin density of 
states of Ni site in each magnetic sublattice. An e^ band 
lies on the energy range Ep — 0.6 eV ~ Ep + 1.2 cV. The 
numbers of states in the energy ranges Ep — 0.6 eV ~ Ep 
and Ep + 0.1 eV ^ Ep + 1.2 eV are 16, respectively, 
which are the one occupied and one vacant Cg states per 
Ni ion. Actually these Cg orbitals extend over surround- 
ing oxygen sites from Ni ions due to strong hybridization 
between Ni Cg and O p orbitals. The extended occupied 
eg orbital has a 60% weight on p orbitals on surround- 
ing six oxygens, a 10% weight on an individual oxygen. 
Therefore, one would establish a model where one oc- 
cupied Bg state with majority spin locates at the top of 
the valence bands, and it hybridizes strongly with the p 
states on nearby O ions. This is the molecular orbital a* 
state. Then one can assign all nickel ions in NdNiOa 
to be trivalent Ni^+ (^296^), even though the Ni ion is 
not truly ionized by -|-3 charge. In the projected density 
of states of Ni ion site, one observes a large amount of 
eg states at the bottom of the d bands , which are the 



bonding states between Ni d and O p, corresponding to 
the a states in the molecular orbital picture, [ p^ 

FIG. 2. Isometric surface of the spin density — pi 
— ±0.01 states/(atomic unit)'' in the energy range 
Ep — 0.6 eV ~ Ep in a thin slab parallel to (001) plane. 
Yellow and blue surfaces indicate the plus and minus spin 
densities, respectively. Black solid lines stand for Ni-O bonds 
and red lines unit cell. 

The eg band in the range Ep — 0.6 eV ^ Ep does not 
show any orbital ordering. In fact, off-diagonal elements 
within the eg subblock of the occupation matrix {n„„„/} 
are zero and the diagonal elements are identical. The ab- 
sence of the orbital ordering may be consistent with the 
fact that the [101] axis has three-fold rotational symme- 
try in the present spin configuration, once the distortion 
is neglected. Due to this pseudo three-fold rotational 
symmetry, the basis orbitals of Eg representation of the 
trigonal group D^d is a good basis set and those derived 
from eg orbitals are Lpu+ — -^{fsz^-r^ + ^fx'^-y^) and 
ipu- — -^{fsz^-r^ — ^fx^-y^)- Therefore, there is no dif- 
ference between the occupancies of ip^^2_^2 and fx^-y^- 

Figures || is the spatial profiles of spin densities p| — pi 
in the energy range Ep — 0.6 eV ^ Ep {a* state). Only 
one spin component of O p orbitals on the FM bond is 
bridging between two Ni Cg orbitals, while both spin com- 
ponents on the AFM bond couple with Ni eg orbitals of 
respective spins. Consequently, oxygen ions on the AFM 
bond have more charge in this energy range than oxygen 
on the FM bond. This is the realization of oxygen-site 
charge-ordered state discussed by Mizokawa et al. in the 
framework of the Hartree-Fock calculation. ||l^ However, 

» FM 

the charge difference between oxygen sites of Ni ' - O -Ni 

AFM 

and Ni' - O -Ni-^ in a* states is mostly compensated by 
hybridized a state at lower energies. Besides, all oxygen 
ions have no local magnetic moment. 

More significant results seen in Fig. ^ may be the p 
orbital ordering on oxygen sites. The magnetic space 
group is P(,a and its unitary part is Pa. The unit cell 
and lattice primitive vectors are not identical to those in 
Pbnm. The spin density in Fig. ^ is consistent with this 
magnetic group P{,a, neither higher nor lower than this. 
The symmetry lowering of the unitary part Pa is the ori- 
gin of the opening the band gap at Ep in the majority 
spin band, despite to the absence of the orbital ordering 
on Ni sites (i.e. symmetry driven band gap). Therefore, 
one can conclude that the origin of the insulator phase 
at low temperatures in NdNiOs is the characteristic spin 
density on the oxygen sites or, equivalently, the orbital 
ordering there. The energy gap Eg = 0.11 eV is due 
to the symmetry of charge order. In fictitious ideal cubic 
structure without distortion or tilting of NiOg octahedra, 
the system becomes metal whose valence and conduction 
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bands touch at points with each other. |T^] The structure 
of <j* bands is insensitive to the value of U . The value 
of the band gap is unchanged down to [/ = 4 eV. This is 
because the gap is driven by the symmetry. 

YNiOs: One should expect larger distortion in YNiOa 
because YNiOa is the typical system with the small toler- 
ance factor. |^,|| There are two different crystallographic 
Ni sites, and the distances from each Ni ion to surround- 
ing O ion are different by from one type of Ni 
ion to the other type. Calculated self-consistent solu- 
tion is that with the apparent charge disproportionation 
and no orbital polarization. Since Ni^+ site has no spin 
magnetic moment, the two spin configurations discussed 
in NdNiOa become identical and the spin configuration 
of YNiOa is uniquely determined. The magnetic space 
group is Pb2i if the spins of Ni^+ ions are non vanishing 
and, on the contrary, Pb2i/a if the spins of Ni^'^ ions are 
zero. The latter symmetry Pf,2i/a is actually the case. 

Figure H shows the projected density of states at Ni ion 
sites. The system is insulator with a gap Eg = 1.03 eV. 
The resultant magnetic moments for half of Ni ions are 
±1.5/iB within the atomic sphere of a radius 2.52 au, 
namely divalent ions Ni^+ (^ig^g); of larger NiOg octa- 
hedron and zero for another half of Ni ions within the 
atomic sphere of a radius 2.51 au, namely quadrivalent 
Ni^+ (t^g); of smaller octahedron. The experimentally 
observed magnetic moments are ±1.4^^ for Ni^+ ions 
and ±0.7/xb for Ni'^^ ions. Q The discrepancy may be 
due to a possible non-coUinear spin order. The number 
of states in the energy range E-p — 0.47 eV ~ E-p eV is 16 
and dominant weight on Ni^+. Because the number of 
Ni^+ in the 2x1x2 cell is 8, these states are assigned as 
two Cg states mainly on Ni^+ and surrounding oxygens. 
The eg orbitals of Ni'^^ is lifted in the higher energy range 
(> E-p -\- 1 eV) without spin polarization. 

A large amount of Cg orbitals in Ni^+ ions locates at the 
bottom of the d bands , in the range Ep — 6.7 eV ^ Ep — 
5.8 eV for Ni^+, and in the range E'p - 5 eV i?F - 4 eV 
for Ni'*+, stabilizes the bonding between Ni'^^ ions and 
oxygen ions. From these facts, one can establish a model 
that deep u molecular orbitals (one per both Ni^+ and 
Ni^+) stabilize the system, and other two Cg states (cr* 
states) per one Ni^+ ion located at Ep — 0.47 eV ~ Ep. 



FIG. 3. Projected density of states of (a) divalent Ni^^ and 
(b) quadrivalent Ni**^ in YNiOa. The energy zeroth is set at 
the top of the occupied states E^. The local coordinate axes 
are chosen in the same manner for those in Fig. 1. 

The charge disproportionation is mainly due to the 
crystal field effect. The low spin state in Ni''"'" or Ni'^"'" 
ion is energetically unstable in small Dq case and the 
ground state multiplet of Ni^"*" is '^A2 (tfgCg) for all arbi- 
trary values of Dq. |l^ Therefore, the small tolerance fac- 
tor causes two different Ni sites, compressed Ni"*"*" (large 



Dq) and dilated Ni^"*" (small Dg), rather than uniformly 
dilated Ni'^"'" ionic states. ||l^ The standard values of 
Dq/B [B is the Racah parameter) for Ni^+ is presum- 
ably around 1.0. Once one estimate the crystal field ef- 
fects from the Ni-0 bond lengths dNi-o, the difference 
of lODg on two Ni ion sites is presumably about 20%. 

Two narrow eg bands can be seen in the energy range 
Ep - 0.47 eV - £;f - 0.23 eV and Ep - 0.23 eV ~ Ep 
in Fig. ||a. Two Ni^"*" Cg states are spatially extending 
over wide area and not only hybridizing with the nearest 
neighbor O ions but also extending over the nearest Ni''^"'" 
ions. This situation is well depicted in the spin density. 
Figure ^ shows the isometric surfaces of the spin den- 
sity in the range of Ep - 0.47 eV ~ iJp - 0.23 eV. The 
d-wavefunctions on Ni^"*" are extending over the (101) 
plane, and antiferromagnetically coupled with other Ni^"*" 
ions. Charge in each Ni ion is compensated in YNiOs as 
in NdNiOa and the difference between total charges in 
the muffin tin spheres on Ni^"*" and Ni''^"'" sites is very 
small, equals to 0.03. This variation is the same order as 
that of oxygen ions in NdNiOa. However, charge dispro- 
portionation is coupled with the lattice distortion, where 
larger oxygen octahedron is surrounding Ni^"*" , and stabi- 
lizes the lattice system in YNiOa. Therefore, diffraction 
experiment can detect the charge ordering in the yttrium 
system easier than in the neodymium system. 



FIG. 4. Isometric surface of spin density p| — pi 
= ±0.002 states/ (atomic unit)^ in the energy range of an Cg 
band {E-p - 0.47 eV~ E-p - 0.23 eV) in a thin slab parallel to 
(101) plane. Yellow surface means the sign of the spin density 
is plus, blue one means minus. 

In conclusion, we have studied two typical antifer- 
romagnetic insulating phase of RNiOa, NdNiOa and 
YNiOa, by using the LSDA+U method. A possibihty 
of the Jahn- Teller distortion is excluded by the group 
theoretical consideration. No orbital order on Ni sites is 
observed in both systems, but two different types of or- 
dering are observed. In small distorting RNiOa such as 
NdNiOa, oxygen sites shows the orbital ordering and this 
is the origin of the gap opening in NdNiOa . In large dis- 
torting RNiOa such as YNiOa , the charge disproportion- 
ation 2Ni'^+ ^Ni^++Ni'*"'" occurs. The charge ordering 
mechanism can explain the stabilization of the asymmet- 
ric alignment of the local magnetic moments around each 
nickel site in both systems. 

We may add the final comment about effects of 
electron-electron correlation in RNiOa • The widths of the 
calculated Cg bands D by the LDA are 2.82 eV(YNiOa), 
2.87 eV(NdNi03), and 3.18 eV(LaNi03). Then, the ra- 
tio of U /D may be estimated as 1 : 0.98 : 0.89, assuming 
the value of the Coulomb repulsion U is common for all, 
and one can see a large reduction of U/D in LaNiOa, 
which may be the key parameter for the difference of 
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the ground states of these perovskite nickelates. NdNiOa 
shows the anomalous M-I transition. Q LaNiOs is pre- 
sumably an anomalous metal of strongly correlated elec- 
trons JlGt and could not be treated within the framework 
of the LSDA+U method. 
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